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1.  Introduction 


Surface  boundary  layer  processes,  in  particular  the  turbulent  kinetic  energy 
budget  and  the  surface  energy  balance,  are  cohesively  linked  to  the  vertical 
profiles  of  wind,  temperature,  and  specific  humidity  that  characterize  the 
mean  flow  and  turbulent  structure  of  approximately  the  first  100  m  of  the 
atmosphere.  The  shapes  of  the  vertical  profiles  are  deemed  to  be  controlled 
by  the  nondimensional  parameters  of  the  dynamic  similarity  theory  of 
Obukhov.  [1]  Priestley  [2]  suggested  that  the  Richardson  [3]  number  was 
the  sole  controlling  parameter  of  similarity  near  the  earth’s  surface.  Two 
significant  entities  that  have  a  profound  e^ect  on  dynamic  similarity  and  the 
surface  energy  balance  must  be  added  to  Priestley’s  corollary,  the 
aerodynamic  roughness  length,  and  the  latent  heat  flux  via  local 
evapotranspiration  rates.  Surface  roughness  length  has  significant  effects 
on  vertical  wind  shear;  larger  shears  are  associated  with  larger  values  of 
surface  roughness.  Large  wind  shears  and  increases  in  the  vertical  latent 
heat  flux  caused  by  increased  evapotranspiration  have  a  stabilizing  effect  on 
the  atmosphere.  The  addition  of  water  vapor  to  the  atmosphere  reduces  the 
sensible  heat  flux,  which  leads  to  the  stabilizing  phenomenon. 

An  orderly  and  comprehensive  examination  of  the  surface  energy  balance 
requires  detailed  analyses  of  the  vertical  profiles  of  wind,  temperature,  and 
humidity  with  respect  to  atmospheric  stability.  Application  of  dynamic 
similarity  principles  to  the  micrometeorological  characteristics  of  the 
surface  boundary  layer  is  interwoven  with  the  nature  of  turbulent  flow  in 
which  the  instantaneous  orthoginal  velocities  exhibit  irregular,  random 
fluctuations  so  that  only  the  statistical  properties  can  be  subjected  to 
analyses.  Thus,  profile  analysis  is  generally  based  on  averages  formed  over 
some  arbitrary  time  interval  and  at  preselected  measurement  levels  on  an 
instrumented  tower  or  mast.  Analytical  arguments  concerning  the  best 
approaches  to  defining  the  profiles  are  undertaken. 
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2.  Buoyancy  in  the  Surface  Layer 


Buoyancy  is  an  atmospheric  property  of  an  air  parcel  that  allows  it  to 
ascend  and  remain  suspended  in  a  compressible  fluid  such  as  the 
atmosphere.  In  the  atmosphere,  a  buoyant  force  on  a  parcel  of  air  may  be 
directly  attributed  to  a  local  increase  in  temperature  resulting  from  the 
heating  of  the  atmosphere  at  the  air-earth  interface  by  conduction,  which 
immediately  leads  to  a  vertical  energy  transfer  by  convection.  The  work 
done  by  buoyancy  against  gravity  results  in  a  decrease  with  height  of 
potential  temperature,  a  characteristic  of  an  unstable  atmosphere. 

Conversely,  during  the  nocturnal  hours,  in  the  absence  of  incoming  solar 
energy,  the  kinetic  energy  of  the  turbulent  fluctuations  will  decrease, 
because  of  the  extraction  of  energy  from  the  mean  motion  by  the  Reynolds 
stresses.  The  atmosphere  becomes  quiescent  and  the  mean  flow  approaches 
laminar  flow  conditions.  The  atmosphere  is  said  to  be  stable. 

The  stability  of  the  atmosphere  may  be  expressed  as  the  ratio  of  mechanical 
production  to  the  convective  production  of  energy.  This  ratio  is  generally 
known  as  the  gradient  form  of  the  Richardson  [3]  number  or  the 
Obukhov  [1]  scaling  length.  Richardson,  in  the  development  of  his 
criterion  of  turbulence,  did  not  frilly  consider  the  effects  of  water  vapor  on 
atmospheric  stability.  Obukhov  also  omitted  this  important  detail  in  his 
arguments  leading  to  the  dynamic  similarity  of  flows  theory. 
Lumley  and  Panofsky  [4]  discussed  the  importance  of  water  vapor  to  the 
dynamics  of  the  atmosphere.  Busch  [5]  redefined  the  Obukhov  length  L 
with  respect  to  the  affect  of  water  vapor  and  the  latent  heat  flux  on 
atmospheric  stability  as 


P9  ^  0^ 

kgiH  +  0.07S££)  kgd* 


(1) 


where 

u.  =  a  scaling  or  friction  velocity 

Cp  =  the  specific  heat  of  air  at  constant  pressure 

p  =  density 
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=  Karman’s  constant 
=  gravitational  acceleration 
=  vertical  sensible  heat  flux 
=  latent  heat  of  vaporization 
=  evaporation  rate 
=  latent  heat  flux 
=  mean  potential  temperature 
=  virtual  potential  temperature 
=  virtual  scaling  temperature. 


By  definition,  the  scaling  length  L  is  related  to  the  gradient  Richardson 
number,  Rj,  by 


z 

L 


(2) 


where 

z/L 

=  Monin  and  Obukhov  [6]  scaling  ratio 

Kh 

=  exchange  coefficient  for  heat 

Km 

=  exchange  coefficient  for  momentum 

0M 

=  dimensionless  lapse  rate. 

Based  on  the  O’KEYPS*  function,  it  may  be  demonstrated  that  in  unstable 
flow  z/L  =  Rj,  because  and  in  a  thermally  stratified  stable 

regime,  z/L  =  Rj  because  =  0m-  Therefore,  in  an  unstable  regime. 


£ 

L 


Zg  aln  Z 


(3) 


*0’KEYPS  is  an  acronym  formulated  by  Panofsky  [7]  and  Yaglom  [8]  and  recognizes  the 
major  contributors  to  the  development  of  contemporary  dynamic  similarity  theory 
Obukhov,  [1]  Kazanski  and  Monin,  [9]  Ellison,  [10]  Yamamoto,  [11]  Panofsky,  [12]  and 
Sellers.  [13] 
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where 


V  =  mean  horizontal  wind  speed 

z,  =  geometric  mean  height  of  the  layer  Zj  -  Zp 

For  stable  mean  flow  conditions,  Hansen  [14]  found  that 
=  (1  +  15  Rj),  allowing  z/L  to  be  expressed  as 

^  =  R,  (1  +  15  R^).  (4) 

Xf 


The  use  of  virtual  potential  temperature  for  surface  boundary  layer  physical 
processes  is  basic  to  proper  and  accurate  solutions.  Virtual  potential 
temperature  is  the  temperature  necessary  for  dry  air  to  produce  the  actual 
density  at  ambient  pressures,  expressed  as 

e,  . - 5 -  (5) 

''  1  -  0.379  elp 


or 

=  6(1  +  0.61  4)  (6) 


where 

e  =  partial  pressure  of  water  vapor 
p  =  atmospheric  pressure 
q  =  mean  specific  humidity. 


7 


3.  The  Mechanics  of  Vertical  Profiles 


Investigation  and  examination  of  surface  boundary  layer  phenomena  usually 
involve  the  collection  of  representative  experimented  data.  In  the  interest 
of  economy,  not  only  in  collecting  experimental  data,  but  also  in  processing 
samples,  vertical  profiles  are  usually  observed.  This  means  data  are 
analyzed  using  bulk  aerodynamic  methods. 

Traditionally,  the  instrumentation  of  towers  and  masts  has  been  at  heights 
above  the  surface  established  by  application  of  the  law  of  natural  growth. 
This  results  in  measurements  at  instrument  levels  spaced  in  a  geometric  or 
logarithmic  progression.  It  also  suggests  that  the  vertical  gradients  are 
tangent  to  the  profiles  at  the  geometric  mean  height  of  each  layer.  Rachele, 
Tunick,  and  Hansen  [15]  investigated  the  validity  of  this  premise  via  the 
mean  value  theorem  of  calculus;  finding  that  the  true  point  of  tangency  was 
not  at  the  geometric  mean  height  of  a  layer,  but  8.2  percent  above 
zg  =  (z,  •  z,)‘^.  Equation  (3)  for  Rj  is  now  written  as 


1 

e  (av)2 


AZ 


(7) 


because,  from  the  mean  value  theorem. 


z  = 


A  z 

A  In  z 


(8) 


This  is  only  of  academic  interest  because  values  of  R^  calculated  for  heights 
of  zg  or  z  translate  into  identical  values  of  L.  The  important  point  is  that 
data  analysis,  using  either  the  law  of  natural  growth  or  the  mean  value 
theorem,  yields  cogent  estimates  of  atmospheric  stability. 

The  stochastic  nature  of  turbulent  flow  implies  that  statistical  methods  must 
be  invoked  to  temporarily  smooth  the  data  for  each  tower  or  mast  level. 
Customarily,  the  Eichardson  numbers  or  Monin-Obukhov  scaling  ratios  are 
formed  at  this  point  in  the  analysis.  An  inspection  of  the  stability 
parameters  usually  reveals  a  nonuniform  distribution  with  height,  and  on 
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conversion  to  Obukhov  lengths,  a  variable  L  with  height.  Dynamic 
similarity  requires  a  constant  scaling  length  throughout  the  surface  layer. 
The  variability  of  L  with  height  can  be  attributed  to  production  and 
dissipation  of  mechanical  and  convective  energy  and  eddies  of  several 
orders  of  magnitude  in  size  impinging  on  or  completely  bypassing  tower 
instrumentation.  Irregular  buoyancy  and  subsidence  effects  on 
anemometers  or  temperature  and  humidity  sensors  yield  nonstationary 
appearing  mean  profile  data. 

Smoothing  the  profile  heights  can  be  accomplished  in  various  ways.  The 
simplest  method  is  smoothing  by  eye  (graphing  the  averaged  data  points  as 
a  function  of  height)  and  hand  fitting  a  smooth  curve  to  the  points.  Values 
picked  from  the  curve  are  used  to  recalculate  the  stability  parameters.  A 
better  fit  can  possibly  be  determined  using  a  technique  attributed  to 
Businger  et  al.  [16]  For  unstable  flow  conditions,  differentiated  second- 
order  polynomials  in  Inz  are  fitted.  For  the  stable  regime,  the  polynomial 
a  +  blnz  =  cz  was  used.  The  smoothed  values  for  each  level  of  interest 
are  used  to  evaluate  the  stability. 

Additional  smoothing  of  the  stability  parameters  is  required  to  provide  best 
estimates  of  the  Obukhov  length.  Lettau  [17]  introduced  the  concept  of  the 
approximate  height  derivative  of  the  gradient  Richardson  number,  that  may 
be  used  as  a  smoothing  technique,  and  is  given  by 


L-‘ 


2(4 


(9) 


where 

i  =  1,  2,  3 

which  is  valid  for  unstable  conditions.  For  the  stable  regime 

L-'  =  (10) 

S(z) 


after  determining  z/L  from  equation  (4).  This  technique  is  referred  to  as 
a  stationarity  forcing  function. 


4.  Dynamic  Similarity  in  the  Surface  Layer 


The  O’KEYPS  designation  has  been  reserved  to  describe  the  functions 
associated  with  the  linear-quartic  profiles  in  unstable  flow.  We  chose  to 
use  the  acronym  as  an  all-encompassing  descriptor  for  both  the  unstable  as 
well  as  the  log-linear  stable  regime.  The  O’KEYPS  originators,  however, 
were  not  the  sole  contributors  to  the  investigation  of  stable  mean  flow 
conditions  in  the  surface  boundary  layer;  there  were  many  others,  including 
Webb,  [18]  Oke,  [19]  and  Hicks.  [20] 

During  this  investigation  concerned  with  water  vapor  effects  on  stability, 
it  soon  became  apparent  that  a  number  of  simplifying  shortcuts  could  be 
used  in  the  data  analyses.  These  abridged  procedures  were  unique  enough 
to  deserve  a  separate  moniker  and  were  dubbed  Mariah.*  O’KEYPS  and 
Mariah  yield  identical  results,  but  Mariah  requires  much  less  computational 
effort. 

The  wind,  temperature,  and  specific  humidity  profiles  in  the  surface 
boundary  layer  of  the  atmosphere,  written  in  differential  form  are 


dz 


u* 

kz 


(11) 


dz 


(12) 


*Mariah  has  no  relationship  to  any  meteorological  phenomenon  or  individual,  but  was 
unabashedly  lifted  from  the  Broadway  Musical  and  motion  picture  "Paint  Your  Wagon ", 
in  which  Uncle  Willie  sang 

way  out  there  they  have  a  name  for  rain  and  wind  and  fire,  the  rain 
is  Tess,  the  fire  is  Joe,  and  they  call  the  wind  Mariah. " 
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and 


dq  ^  qi 
dz  kz 


4>w 


(13) 


where 

T*  =  a  scaling  temperature  associated  with  potential  temperature 
q*  =  a  scaling  .tumidity. 

Within  the  framework  of  dynamic  similarity,  it  is  presumed  that  the  transfer 
of  heat  and  mass  are  identical,  so  0*  =  <f>H-  From  the  O’KEYPS  function, 
it  is  easily  demonstrated  that  for  unstable  conditions 

(|)„  =  (1-  15  z/t)V*  (14) 


and 


♦w  =  (1-  15 


Integrating  equations  (11),  (12),  and  (13)  using  equations  (14)  and 
(15)  yields 


K  =  -y  [In  z/z^  +  (z/L)]  (16) 


0  -  00  =  y  [In  z/z^  +  i|fj^z/I)] 


(17) 
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and 


q  -  <lo  =  ^  ^ 


(18) 


The  diabatic  influence  functions  are  given  by 


and 


W)  =  -  /  -p-  d  (zlD 

r*  1  -  <!>*/ 

(^^)  =  "  /  ^ 


(19) 


(20) 


and  solved  using  Benoit’s  [21]  method 


=  - 


1 1n 


(S^  1)  Stp  H- 

(3?  -  1)  &o  ■>  D* 


+  2  [ran’*  (Sf)  -  Tan-'  (2^] 


(21) 


where 


X 

X 


1 


1/4 
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and 


CUL)  =  2  In 


(22) 


where 


1 


15 


L 


The  scaling  temperature  for  unstable  flow  may  be  found  from 


r  = 


e«; 

k  zg 


(23) 


and  according  to  Myrup  [22],  q^,  the  specific  humidity  at  Zq,  is  determined 
from 


9o  = 


1000 


» 

(t.) 

2 

3.74  +  2.64 

V 

iloj 

(24) 


where 

To  =  the  temperature  at  ZIq  in  degree  Celsius. 
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Equations  (16),  (18),  and  (23)  now  can  be  evaluated  for  u*,  q*,  and  T*,  the 
initial  scaling  parameters  required  to  solve  the  surface  energy  balance.  This 
approach  is  laborious  because  it  is  necessary  to  first  evaluate  the  diabatic 
influence  functions  in  equations  (21)  and  (22). 

A  second  approach  is  to  guess  a  value  for  the  Obukhov  length  L  (negative 
for  unstable  conditions  and  positive  for  stable).  Using  the  value  for JL,  a 
solution  is  found  for  equations  (21)  and  (22).  For  known  values  of  V,  6, 
and  q,  u.,  T*,  and  q*  are  evaluated  using  equations  (16),  (17),  and  (18). 
The  values  are  used  to  determine  H  =  -  Cppu.T*  and  ifE  =  -  ifpu*q*  that 
are  substituted  in  equation  (1)  to  obtain  a  new  estimate  of  L.  This  value  of 
L  is  used  to  form  new  estimates  of  u*,  T*,  and  q*.  The  process  is  repeated 
until  the  scaling  constants  converge.  The  second  approach  is  not  very 
efficient,  even  with  machine  processing. 

In  a  thermally  stratified  stable  regime  in  which  z/L  =  Rj 
the  integration  of  equations  (11),  (12),  and  (13)  yields 


V 


u* 

k 


(25) 


0  -  00 


II 

k 


In  ?/Zo  +  P  2 


(26) 


and 


In  ?/Zo  +  P  Y 


(27) 


where 

/3  =  a  variable  according  to  Hansen.  [14] 
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If  ^  =  1  +  j8  z/L,  then  from  equation  (14) 


(28) 


where 

=  (1  +  15  Ri). 

The  scaling  parameters  may  then  be  found  from 


u* 


itAV 
A  Inz  +  p 


(29) 


r  = 


k  aQ 
A  In  z  +  p 


(30) 


and 


9*  = 


k  Aq 


Alnz  +  p  ^ 


(31) 


As  with  the  unstable  case,  a  value  of  L  is  needed  to  obtain  values  of  u*,  T, 
and  q*  using  equations  (29),  (30),  and  (31).  This  is  done  by  employing  an 
iterative  approach  similar  to  the  unstable  case. 


5.  The  Mariah  Approach 


The  Mariah  approach  evolved  from  a  study  by  Hansen.  [23]  From  the 
definition  of  the  Obukhov  [1]  scaling  length  of  equation  (1),  repeated  here 
as 


(32) 


where 


0*  =  r  +  0.61  0^*, 


(33) 


we  may  substitute  equations  (11),  (12),  (13),  and  (32)  in 

equation  (1),  yielding 


L  = 


e, 

g  A  In  ^  (^  0  ()*61  0  ^4) 


(34) 


which  is  valid  in  both  the  stable  and  unstable  regimes,  and  allows  a  direct 
calculation  of  L  from  the  observed  gradients.  However,  in  stable 

conditions  0H  =  0m>  therefore,  the  denominator  contains  0m-  Manipulation 

of  equation  (33)  reveals  that  in  stable  flow 


L  = 


e,  (aF)' 


g  A  /nz  (a  0  +  0.61  0  Ag) 


-  15 


A  Z 

Inz’ 


(35) 
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which  leads  to 


L  = 


60  Z  +  2B 


(36) 


with  z  =  az/a  Inz,  and 


B  = 


2  e,  (aIQ^ 

^  him  (a  6  0.61  6  A^) 


(37) 


a  solution  for  the  Obukhov  length  independent  of  the  dimensionless  wind 
shear. 

The  scaling  constants  u*,  T*,  and  q*  are  calculated  using  the  discrete  forms 
of  equations  (11),  (12),  and  (13), 


and 


u* 


r 


9* 


fa:  aK  _ 

k  hV 

4>j,  hlnZ 

feA0  _ 

k  hV 

aZ  hliiZ 


kz  hjq  _ 

k  tq 

01,  AlnZ’ 

(38) 


(39) 


(40) 


and  while  evaluating  0^  using  equation  (14),  noting  that  z  =  Az/Alnz. 


6.  Discussion 


The  Mariah  approach  is  simpler  than  previous  solutions  involving 
similarity.  Aside  from  directly  estimating  the  Obukhov  scaling  length  L, 
Mariah  produces  results  that  do  not  differ  numerically  from  the  O’KEYPS 
equations.  Mariah  operates  using  measured  gradients,  and  O’KEYPS 
operates  using  chosen  trial  values  and  iterating  until  convergence.  Hand 
calculation  using  the  O’KEYPS  equations  is  tedious. 

Figure  1  shows  a  comparison  of  scaling  specific  humidities  determined 
using  both  schemes.  The  significant  feature  is  the  evaluation  of  q*  in  the 
stable  regime.  During  the  nocturnal  hours,  evaporation  continues  to  the 
point  before  the  humidity  profiles  go  isohumic  and  invert.  Inversion  of  the 
specific  humidity  profiles  is  the  first  step  toward  dew  fall.  The  transfer  of 
water  vapor  from  the  atmosphere  to  the  surface  appears  to  have  little  or  no 
effect  on  the  application  of  dynamic  similarity  to  the  surface  boundary 
layer. 

Figure  2  shows  specific  humidity  effects  on  atmospheric  stability.  An 
increase  of  water  vapor  in  the  atmosphere  in  conjunction  with  larger  latent 
heat  fluxes  reduces  absolute  values  of  the  Richardson  number  or  Obukhov 
scaling  length  4  to  10  percent. 

Brunner  [24]  suggests  that  in  the  Planetary  Boundary  Layer  relative 
humidity  can  be  considered  a  constant  in  approximately  ^e  first  100  m  of 
the  atmosphere.  Thus,  relative  humidities  calculated  from  wet  and  dry  bulb 
temperatures  at  shelter  heights  can  be  easily  translated  to  vapor  pressures 
and  specific  humidities  at  tower  heights.  Although  it  is  not  the  most 
exacting  method,  it  will  provide  the  moisture  corrections  needed  for  surface 
energy  balance  studies. 
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Figure  1.  Comparison  of  scaling  specific  humidities  determined  using  both  schemes. 
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7.  Conclusions 


The  Mariah  approach  to  similarity  theory  profile  analysis  is  simpler  than 
that  of  the  O’KEYPS  methodology,  especiidly  if  the  more  formal  analyses 
involving  the  diabatic  influence  functions  (z/L)  and  used 

to  evaluate  the  scaling  parameters  u*,  T*,  and  q*  in  the  unstable  regime. 

O’KEYPS  and  Mariah  appear  to  function  well  in  the  stable  regime  — 
especially  after  inclusion  of  water  vapor  in  the  models.  The  approach  to 
an  isohumic  state  and  the  inversion  of  the  specific  humidity  profiles  seems 
to  have  little  or  no  effect  on  the  profile  analyses  schemes,  indicating  that 
expansion  of  the  Mariah  scheme  by  estimating  specific  humidity  profiles  for 
shelter  height  relative  humidities  is  highly  feasible  for  the  entire  stability 
spectrum. 
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